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We have carried out a computational study of the inelastic electron tunneling spectrum (IETS) of 
the two vibrational modes of a single hydrogen atom on a Cu(100) surface in a scanning tunneling 
microscopy (STM) junction. This study addresses key issues about vibrational assignment and line 
shape of observed peaks in IETS within the framework of density functional theory calculations and 
the Lorente-Persson theory for STM-IETS. We argue that the observation of only a single, broad 
peak in the STM-IETS [L. J. Lauhon and W. Ho, Phys. Rev. Lett. 85, 4566 (2000)] is not caused by 
any symmetry restrictions or any cancellation between inelastic and elastic vibrational contributions 
for one of the two modes but is due to strongly overlapping superposition of the contributions from 
the two modes caused by the rather large instrumental broadening and the narrow vibrational energy 
separation between the modes. In particular, we find that this broadening and the large asymmetry 
of the vibrational line shapes gives rise to substantial apparent vibrational energy shifts of the two 
modes and decrease their apparent energy separation. 
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A most important advance in surface and nano science 
has been the realization of single molecule vibrational 
spectroscopy and microscopy at metal surfaces using a 
scanning tunneling microscope (STM)i. The unique ca- 
pabilities of the STM are based on two key ingredients: 
(1) the highly localized tunneling of electrons in the space 
between the tip and the sample, resulting in atomic scale 
resolution in imaging; (2) the non-adiabatic coupling be- 
tween tunneling electrons and the vibrations in the STM 
junction giving rise to characteristic signatures in the 
tunneling current I at biases V corresponding to vibra- 
tional energies. Usually, these thresholds are identified 
as peaks in the jkrziV) signal. The observed vibrational 
energies provide unique information about the chemical 
identity of single adsorbed species and can also be used 
to discriminate among various isotopes. To fully exploit 
this new vibrational spectroscopy, one needs physical in- 
sights from theory about the strengths, line shapes, and 
lateral spatial distributions of the vibrational signatures 
in the j^k(V) signal. 

The theory of inelastic electron tunneling spectroscopy 
(IETS) from vibrations in a STM junction has been 
developed at various levels of sophistication by sev- 
eral groups 2 * 3 *^. One important advance was made by 
Lorente and Persson^, who generalized the widely used 
Tersoff-Hamann theory for elastic tunneling to inelastic 
tunneling from vibrations and implemented the theory in 
density functional theory (DFT) calculations. As demon- 
strated by a direct comparison of the results from the 
LP theory with IETS experiments, this theory is able to 
reproduce the observed strengths of various vibrational 
modes for single molecules on a copper surface^i^. In 
particular, the LP theory suggested that many modes 
were not observed due to a cancellation between inelastic 
and elastic contributions to the tunneling current. Fur- 
thermore, the results of the LP theory suggested a sym- 



metry selection or rather propensity rule for IETS relat- 
ing the symmetry of the vibrational mode, the spatial 
distribution of the inelastic signal and the tip orbital 9 . 
The existence of such a rule is of key importance in the 
assignment of various peaks in the IETS to vibrational 
modes. 

A most simple system that raises issues about the role 
of symmetry and the assignment in IETS is provided by H 
and D atoms adsorbed on Cu(100). In these experiments 
by Lauhon and He 10 , IETS experiments were used in a 
beautiful manner to discriminate between the H and D 
atoms, enabling the study of the isotope effect on surface 
diffusion. However, they observed only a single broad 
peak that was attributed to the perpendicular vibrational 
mode. This observation raises questions about the origin 
of the large broadening of the observed peak and why the 
parallel mode is not observed. 

In this letter, we present a theoretical study of the 
IETS of H and D atoms on a Cu(100) surface based on the 
LP theory and DFT calculations. In particular, we clarify 
the origin of the observed broad single peak in the IETS 
in terms of strongly overlapping contributions from both 
the perpendicular and parallel vibrational modes, caused 
by the extrinsic broadening and the intrinsic asymme- 
try of the vibrational line shape. Before presenting and 
discussing our results we begin by a short review of the 
necessary concepts and ingredients of the theory and the 
calculations of IETS. 

The LP theory of STM-IETS is based on the simple 
and physically transparent Tersoff-Hamann theory for 
imaging by elastic electron tunneling 11 . This theory is 
based on a simple tip model, in which an emitted elec- 
tron is approximated by an s wave, and the Bardeen ap- 
proximation for tunneling. Using these approximations 
one obtains that the differential conductance, 7^7 (U), at 
small biases, V, and low temperatures is simply deter- 



mined by the local density of sample states (LDOS) at 
the position tq of the tip apex as, 



— (F) ocp(ro,e F + e^), 



(1) 



where tp is the sample Fermi energy. In an one-electron 
approximation, as provided by density functional theory 
p(r, e) is given by 



P(r, e) = ^2\(ro\ ip a ) | 2 5(e a - e) 



(2) 



where (fo | i> a ) is an one-electron wave function with 
energy e a . 

The LP theory is based on the observation that the 
TH theory for elastic tunnelling can be directly gener- 
alized to inelastic tunnelling from adsorbate vibrations 
by considering the many-body LDOS for the electrons 
interacting with the adsorbate vibrations. This result 
is based on a few physical assumptions that are fulfilled 
in most cases. First, the electron- vibration coupling is 
assumed to be short-ranged and limited to the sample. 
Second, the vibrational relaxation rate should be much 
more rapid than the tunnelling rate, so that the vibration 
is nearly equilibrated. Third, the electron-vibration cou- 
pling is weak so that the vibration-induced LDOS, Ap 
can be evaluated by first order perturbation theory. This 
non-adiabatic coupling gives rise to two distinct contribu- 
tions to Ap — one from opening up an inelastic channel 
for tunnelling and one from its influence on the elastic 
tunneling channel. 

The channel for inelastic tunnelling from a vibration 
with vibrational frequency Q opens up when the bias V is 
larger than hQ/e and increases in general the tunnelling 
current. This threshold in / results in a positive peak in 
the line shape function 



L (V) 
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where pbg is the LDOS associated with the background 
differential conductance (dI/dV)b s . The integrated 
strength rii ne i(ro) is given by the inelastic fraction of the 
tunnelling electrons with an energy larger than hfl. In 
the quasi-static limit, O — > 0, this fraction is given by the 
vibration-induced deformation of the vacuum tails of the 
onc-clcctron wave functions as 
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(4) 



Here (ro | 6ip a ) is the change in (fo | ip a ) by an rms 
displacement ^jTi/2m£l of the vibrational mode. 

The non-adiabatic coupling of the tunnelling electrons 
with the adsorbate vibration affects also the elastic tun- 
neling channel at the threshold for the inelastic tunneling 
channel through the Pauli exclusion principle. This will 
show up as a decreased integrated strength and an asym- 
metry of the vibrational line shape, as given by, 



Lo{V) 



r 1 (r )( 1 /2) + X (r )(eV-hn) 
Tr((eV-hny + (j/2) 2 ) 



(5) 



Here the total integrated strength r](f ) = r]i ne i(fo) + 
7y c i(ro) has a negative contribution from the elastic chan- 
nel given by, 

%i(r ) = -—V |Im(W a | fo}) | 2 S(e a - e F ). (6) 
The asymmetry parameter x( r o) is given by, 



X(r ) = Y]Re((f o |^a))Im(((5'0Q I r ))6(e a - ep) 

(7) 
The remaining parameter determining the line shape in 
Eq. JSJ is the vibrational relaxation rate 7. 

For high frequency adsorbate vibrations on metal sur- 
faces such as vibrations of H (and D) atoms on a Cu sur- 
face, 7 is dominated by the absorption and emission of 
single electron-hole pairs. The electron-hole pair contri- 
bution 7 c h to 7 can be calculated also in the quasi-static 
limit from, 12 



7eh 



iirfl J2\(^\ Sv I V/3) | 2 S(e a - e F )S{e - e F ) 
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Here Sv is change of the one-electron potential induced 
by an rms displacement of the vibrational coordinate. 

The density functional theory calculations of the elec- 
tronic and geometric structure of an H atom on a Cu(100) 
surface were carried out using Vienna ab-initio simulation 
package VASP 1 ^. The electron ion core interactions were 
handled by the plane wave projector augmented wave 
(PAW) method^ and exchange-correlation effects by a 
generalized gradient approximation^. Because of the 
large lateral extension of the STM image and IET sig- 
nals, we have to use a large 6x6 supercell with Cu atoms 
in six layers, separated by a vacuum region corresponding 
to five layers of atoms. For this large supercell it was suf- 
ficient to sample the Brillouin zone with a 2x2x1 mesh. A 
cut-off energy of 20 Ry for the plane waves was found to 
be sufficient. In the calculations, the H/Cu(100) system 
was first geometrically optimized until forces on each ion 
were smaller than 0.02 eV/A. The technical details about 
the calculations of the electron- vibration matrix elements 
and the vacuum tails of the wave functions, and the han- 
dling of the discrete set of electron states can be found 
in Ref. 6 . 

The experimentally determined four-fold hollow ad- 
sorption siteifi was chosen as the initial configuration in 
the calculations. This site has been shown to be the ener- 
getically preferred site also in previous DFT studiesi£*ii 
After full structural relaxation, the vertical equilibrium 
distance of hydrogen from the uppermost copper layer 
was 0.58 A. The nearest neighboring copper atoms to 
the hydrogen atom relaxed 0.04 A outwards from the 
surface layer and less than 0.02 A parallel the surface 
away from the hydrogen adsorption site. Relaxations of 
the other copper atoms were found to be much smaller. 
The calculated adsorption energy -0.182 eV relative to 
the binding energy of the H atom in the free H2 molecule 
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FIG. 1: Topographical image of constant LDOS of an H atom 
adsorbed in the hollow site on a Cu(100) surface. Az is the 
change in tip-surface distance from its value of 6.26 A away 
from the H atom. The positions of the Cu atoms in surface 
layer are indicated by the open diamonds. 



(E ads = £h/Cu(ioo) - -Ecu(ioo) - !-Eh 2 ) is in good agree- 
ment with the value -0.179 eV obtained by previous DFT 
calculations^. 

In Fig. Q] we show the topographical image of con- 
stant LDOS of the H adatom at a typical value for the 
tip-surface distance. The main characteristic feature of 
this weakly corrugated image is the 0.05A depression at 
the adsorption site surrounded by four 0.06A protrusions 
centered close to the nearest neighbouring hollow sites in 
the [001] and [010] directions. The shape of the image 
with the single depression surrounded by the four pro- 
trusions is robust for all tip-surface distances between 5 
and 10A. However, the depression at the hydrogen site 
decreases in depth with increasing tip-surface distance 
zq: the depth decreases from 0.07A to 0.03A when z$ 
increases from 5 to 8A. Note that the overall size of the 
affected area, about lOxlOA 2 , is very large compared to 
size of hydrogen atom. All the features in the simulated 
STM image are in excellent agreement with the experi- 
mental imageiS. 

The origin of the protrusions in the STM image is 
not simply due to the outward relaxation of the near- 
est neighboring copper atoms or any H-induced state at 
the Fermi level. The protrusions still exist in the sim- 
ulated image of the H atom on the unrelaxed, bare Cu 
surface. In accordance with the common picture of the H 
chemisorption 1 — , the H atom induces an s-like state just 
below the onset of the Cu d band about 7 eV below the 
Fermi level. The small depression in the simulated image 
and the surrounding protrusions then originates from the 
orthogonalization of the metal states at the Fermi level 
against the H-induced state and also from the oscillatory 
screening response by the metal states of the H - entity. 

The H adatom has two vibrational modes which in 
principle can contribute to the IET vibrational spec- 
tra. In the harmonic approximation, the calculated vi- 



brational energies of the perpendicular mode and the 
twofold degenerate parallel mode are 84 and 68 meV, 
respectiveljii^. These values are somewhat different from 
the corresponding values of 76 and 70 xaeSf— (71 and 62 
meV— ) found in earlier density functional calculations 
which take into account anharmonic effects. In these cal- 
culations the vibrational excitation energies were deter- 
mined from the energies of the ground state and the first 
excited states of A\ and E symmetry character in the 
full potential energy surface. These results are in better 
agreement with the experimental value of 70 meV— for 
the dipole active (perpendicular) mode and will be used 
in the detailed comparison with the IET experiments. 

In Fig. El we depict the two-dimensional contour plots 
of the calculated IET line shape parameters rj c \(ro), 
Vmci(ro), r)tot(r Q ) and x(ra), as defined in Eqs. |U Eland 
|7| around the H atom for the two vibrational modes^i. 
For the doubly degenerate parallel modes it is sufficient to 
map out the IET line shape parameters of the mode with 
a vibrational displacement in [011] direction since for the 
other degenerate mode (with a displacement in the [Oil] 
direction) the parameters are simply obtained by rotat- 
ing the maps 90 degrees. The parallel mode parameters 
are then obtained simply as a sum of the contributions 
from the two degenerate modes. 

The vibrational inelastic and elastic signals behave dif- 
ferently for the two vibrational modes. In the case of the 
perpendicular mode, both 7701(^0) and 77inei(^b) are peaked 
on top of the H adsorption site, and cancel each other to a 
large extent. The extreme values for 77 G i(ro) and rji ne i(fo) 
are -1.1% and 2.0%, respectively so that rj to t(fo) has a 
maximum value of 0.9% on top of the H atomSS. In con- 
trast, there is no near cancellation of 7y i(ro) and 77^1(^0) 
for the parallel mode. Thus, the extreme value of 1.4% 
for r/totifo) is stronger than for the perpendicular mode 
despite that the extreme values -0.4% and 1.4% for J7 e i(r ) 
and if me \(fo), respectively, are both weaker than the cor- 
responding values for the perpendicular mode22. Fur- 
thermore, the extreme value of rjtot(ro) from the twofold 
degenerate parallel modes will be twice as large since this 
value is attained at the symmetrical site on top of the H 
atom. Note that the IET line shape parameters depend 
somewhat on the tip-surface distance; maps in Fig.|2]cor- 
respond to a typical distance of 7A2£. 

Although the perpendicular and parallel vibrational 
modes have different symmetry characters (Ai and E), 
the proposed symmetry selection rule in Rcf. 9 or rather 
the symmetry propensity rule for IET does not give any 
definite predictions about the strengths and spatial be- 
haviors of the inelastic and elastic signals for these two 
modes in this case. The effect of the H atom on the states 
around the Fermi level is minute as demonstrated by the 
small corrugation of the calculated topographical image 
of the H atom. This result shows that the LDOS at the 
Fermi level and around the H atom will not be dominated 
by states with a definite symmetry character. Thus there 
will be no significant symmetry restrictions for the IET 
in this case. As shown in Fig. [21 this fact is reflected 
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FIG. 2: The first three rows of the columns show the spatial maps of the calculated elastic, inelastic and total signal strengths 
for the perpendicular (left column) and the parallel mode vibrating in [Oil] direction (right column), respectively. The fourth 
row of figures show the spatial map of the calculated asymmetry parameter for each mode. The H atom is located at the 
center of the images, and the position of the Cu atoms in the surface layer are indicated by open diamonds. The x- and j/-axis 
correspond to the [Oil] and [Oil] crystal directions, respectively. Each image has the same size of 15.4Axl5.4A. The tip-surface 
distance is 7 A in all figures. 



by the spatial maps of rjtot being qualitatively same for 
the perpendicular and parallel vibrational modes. Hence 
these maps cannot simply be used to discriminate be- 
tween these two modes. 

One qualitative difference in the IET from the two vi- 
brational modes is the sign of the asymmetry parameter 
x(?o), which has important implications for the vibra- 
tional line shapes in the IETS. At the H site, x(^b) f° r 
the perpendicular mode is negative with an extreme value 
of -1.2% whereas for the parallel mode it is positive with 
an extreme value of 0.7%. Both these values are com- 
parable in magnitude to rjtot^o), and should give rise to 
vibrational IETS line shapes with characteristic asymme- 
tries for the two modes. To show that, we need first to 
include the calculated vibrational line widths, obtained 
from Eq. [8] into Eq. [S] describing the vibrational line 
shapes. The electron- hole pair line widths j e h are 1.8 
and 1.1 meV for the parallel and perpendicular modes of 
the H adatom, respectively^*. The corresponding relax- 
ation rates 2.7- 10 12 s _1 and 1.6 • 10 12 s _1 are much larger 
than the tunnelling rate I/e ~ 10 9 s _1 corresponding to 
the 0.1 nA current used in the experiments^. Thus the 
current-induced non-equilibrium population of the mode 
should be small justifying one of the assumptions behind 
the LP theory that the vibration should be close to equi- 
librium. 

In Fig. |3 we show calculated vibrational IET line 
shapes for the two modes of the H adatom with the tip 
on top of the H atom. The intrinsic line shapes Lq(V) 
(dashed lines) calculated from Eq. [3] are clearly narrow 
and asymmetric as shown in Figs. EH&) and (b). Thus 
we would expect that the IETS should have two clearly 
resolved peaks in apparent conflict with the single broad 
peak in the IETS observed by Lauhon and Ho. How- 
ever, the theory can be reconciled with experiments by 
taking into account the large extrinsic broadening pro- 
vided by the temperature smearing of the tip and sam- 
ple Fermi distributions, and the modulation voltage2*. 
In the experiments by Lauhon and Ho, the modulation 
voltage had an rms value of 7 mV and the temperature 
was 9 K. Including these broadenings for the two modes 
results in the extrinsic line shape shown as solid lines in 
Figs. EHa)-(c). The large extrinsic broadening and the 
opposite asymmetries of the two line shapes result in 
broad peaks that are significantly red and blue shifted 
for the perpendicular and parallel modes, respectively. 
In addition to decreasing the apparent energy separation 
between the peaks, the broadening is so large that the 
two peaks overlap strongly and can no longer be resolved 
as two peaks. This effect is even more pronounced when 
using the more precise values of 71 meV and 76 meV for 
the vibrational excitations obtained from Ref. 16 . Finally, 
the calculated total strength of the IET from the two 
modes is about 2-3% depending on the distance, which 
is in reasonable agreement with estimate of about 1.1% 
for the observed peak in the IETS. The theoretical iso- 
tope dependence m _1 of rjtot with the atomic mass m is 
consistent with the observed isotope dependence of the 
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FIG. 3: Calculated line shapes of the normalized second order 
differential conductance for (a) the perpendicular mode, (b) 
the doubly degenerate parallel mode, and (c) the sum of the 
modes. The intrinsic line shapes are shown as dashed lines, 
while the line shapes that include temperature and modula- 
tion effects are shown as solid lines. All these line shapes are 
obtained using the calculated vibrational energies of 68 meV 
and 84 meV for perpendicular and parallel modes, respec- 
tively, as obtained in the harmonic approximation. The cal- 
culated line shape based on calculated vibrational energies 
(71 meV and 76 meV for the perpendicular and parallel, re- 
spectively) from Rcf. 16 that takes into account anharmonicity 
and includes temperature and modulation effects is shown as a 
dash-dotted line in (c). The conductances correspond to a tip- 
surface distance of 7 A with the tip on top of the H atom and 
T) = 1.4(0.9)%, x = 0.7(-1.2)%, and 7 = 1.79(1.06) meV for 
the twofold degenerate parallel mode (perpendicular mode), 
respectively. 



intensity for the H and the D atom. 

In summary, we have carried out a computational 
study of the inelastic electron tunneling spectrum (IETS) 
of the two vibrational modes of a single hydrogen atom 
on a Cu(100) surface in a scanning tunneling microscopy 
(STM) junction. This study addresses key issues about 
vibrational assignment and line shape of observed peaks 
in IETS and is based on density functional theory calcu- 
lations and the Lorente-Persson theory for STM-IETS. 
We argue that the observation of only a single, broad 
peak is not caused by any symmetry restrictions or any 
cancellation between inelastic and elastic vibrational con- 
tributions for one of the two modes but is rather a super- 
position of the contributions from the two modes caused 



by the large instrumental broadening and their narrow vi- 
brational energy separation. In particular, we find that 
this broadening and the large asymmetry of the vibra- 
tional line shapes gives rise to substantial apparent vibra- 
tional energy shifts between the two modes and decrease 
their apparent energy separation. 
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